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ABSTRACT 

We present the average abundances of the intermediate elements obtained by performing a stacked analysis 
of all the galaxy clusters in the archive of the X-ray telescope ASCA. We determine the abundances of Fe, Si, 
S, and Ni as a function of cluster temperature (mass) from 1 - lOkeV, and place strong upper limits on the 
abundances of Ca and Ar. In general. Si and Ni are overabundant with respect to Fe, while Ar and Ca are very 
underabundant. The discrepancy between the abundances of Si, S, Ar, and Ca indicate that the a-elements 
do not behave homogeneously as a single group. We show that the abundances of the most well-determined 
elements Fe, Si, and S in conjunction with recent theoretical supernovae yields do not give a consistent solution 
for the fraction of material produced by Type la and Type II supernovae at any temperature or mass. The general 
trend is for higher temperature clusters to have more of their metals produced in Type II supernovae than in 
Type las. The inconsistency of our results with abundances in the Milky Way indicate that spiral galaxies are 
not the dominant metal contributors to the intracluster medium (ICM). The pattern of elemental abundances 
requires an additional source of metals beyond standard SN la and SN II enrichment. The properties of this 
new source are well matched to those of Type II supernovae with very massive, metal-poor progenitor stars. 
These results are consistent with a significant fraction of the ICM metals produced by an early generation of 
population III stars. 

Subject headings: galaxies: abundances — intergalactic medium — supernovae: general — X-rays: galaxies: 
clusters 



1. INTRODUCTION 

Galaxy clusters provide an excellent environment for de- 
termining the relative abundances of the elements. Because 
clusters are the largest potential wells known, they retain 
all the enriched material produced by the member galaxies. 
This behavior is in stark contrast to our own Milky Way 
jTimmes. W ooslev. & Weaver 1995) and many other individ- 
ual galaxies (Henrv & Worthev 1999). The accumulation of 
enriched material in clusters can be used as a probe to study 
the star formation history of the universe, the mechanisms that 
eject the elements into the ICM, the relative importance of dif- 
ferent classes of supernovae, and ultimately the source of the 
metals in the intra-cluster medium (ICM). 

The dominant baryonic component in clusters is the hot gas 
in the ICM, with 5-10 times as much mass as resides in the 
stellar component. The physics describing the dominant emis- 
sion mechanism of the ICM gas is relatively simple. The ICM 
is optically thin, well modeled by a sphere of hydrostatic gas 
in thermal equilibrium, and the high temperatures and mod- 
erate densities minimize the importance of dust. Extinction, 
ionization, non-equilibrium and optical depth effects are min- 
imal. As a result, cluster abundance determinations are more 
physically robust and reliable than those in, e.g., stellar sys- 
tems, H II regions, and planetary nebulae. The hot gas emits 
dominantly by thermal bremsstrahlung in the X-ray band, and 
the strong transitions to the n=l level (K-shell) and to the n=2 
level (L-shell) of the H-like and He-like ions of the elements 
from carbon to nickel also lie in the X-ray band. This makes 
the X-ray band an attractive place for elemental abundance 
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determinations. 

Early X-ray observati ons of galaxy clusters 
(Mitchell. Culhane. Davison. & IveslfT97 6: Serlemitso s et al.l 
1977) showed that the strong H- and He-like iron lines at 6.9 
and 6.7 keV could lead to a value for the meta l abundance 
in clu sters. Later results ( Mushotz kv et a 
Il983h derived from iron line observations showed that 
clusters had metal abundances of about 1/3 the solar value. 

The improved spectral resolution and large collecting area 
of the ASCA X-ray telescope brought new power to studies 
of cluster metal abundances. In particular, the improved 0.5- 
lO.OkeV energy range of ASCA allowed for better spectro- 
scopic fits to clusters than was possible with ROSAT, which 
had an upper limit of 2.5 keV. Mus hotzkv et alJ 0996) stud- 
ied four bright clusters at temperatures such that strong line 
emission is present, and provided the first measurements of 
elemental abundances oth er than iron since the initial results 
from Einstein jM ushotz kv et all (1 981: Becker etal. 197§ 
'Rothe nflug. Vigroux. Mushotzkv. & Holt 1984). Their mea- 
surements of silicon, neon and sulfur were interpreted as high 
abundances of the a-elements in clusters. This result sug- 
gested that type II supernovae (which produce much higher 
a element yields than SN la) from massive stars are responsi- 
ble for a significant fraction of the metals in the ICM. (Type 
la supernovae (SN la) produce high yields of elements in the 
iron peak, while Type II supernovae (SNII) produce yields 
rich in the a elements Si, S, Ne, and Mg.) Later work by 
iFukaza wa ( 1997) showed that clusters are more metal en- 
riched in their centers, and that the Si/Fe ratio is about 1.5- 
2.0 with respect to the solar value. Fukazawa et al. (199^ 
showed how the silicon abundance was higher in hotter clus- 
ters, and confirmed the importance of SN II i n cluster enrich - 
ment. More recently, Finoguenov. David. & Pgnmgi]|j|2222L, 
iFino guenov. Arnaud. & David (i200 1). and Finog uenov et alJ 
(2002) used ASCA and XMM data to show that type la prod- 
ucts dominate in the centers of certain clusters and how type 
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II products are more evenly distributed. The observation by 
IXrnaud et al. ( 1992) that the metal mass in clusters is corre- 
lated with the optical light from early type galaxies and not 
from spirals is also important in determining the origins of 
metals in the ICM. 

In this paper we use the ASCA satellite 
jTanaka. Inoue. & Hoi Ll994) to further constrain the 
abundances of the intermediate elements. Previous authors 
referred to the elements Ne, Mg, Si, S, Ca and Ar as 
a-elements in order to emphasize their supposed similar 
formation mechanism; we will refer to the elements ob- 
servable with X-ray spectroscopy in the ASCA band as 
intermediate elements. This label includes nickel in the group 
and is preferred since the observations will show that the 
a-elements do not act homogeneously as a single class. 

The large database of over 300 cluster observations makes 
the ASCA satellite well suited for a overall analysis of the in- 
termediate element abundances in galaxy clusters. While it is 
not possible to obtain accurate abundances of these elements 
for more than a few individual clusters, we jointly analyze 
many clusters at a time in several "stacks" in order to obtain 
the signal necessary for obtaining the abundances. The rel- 
atively large field of view of the ASCA telescope allows for 
spectroscopic analysis of the entire spatial extent of all but 
the closest clusters, and the moderate spectral resolution en- 
ables abundance determinations from the K-shell and L-shell 
lines. 

iFukazawal Jl997h and other observations of clusters ob- 
tained with Chandra and XMM have show that abundance gra- 
dients are common across the spatial extent of clusters, often 
with enhanced iron abundances in the cluster centers. These 
observations shed valuable light on the source of the metals 
and help discriminate among the mechanisms that enrich the 
ICM. [DeGrandi (2003)^ has shown with BeppoSax measure- 
ments that the centers of clusters (within a radius where the 
density is 3500 times the critical density) have iron abun- 
dances that are enhanced by about 10-20%. These results 
indicate the importance of a physical mechanism in the very 
center of clusters that causes an increase in the central metal- 
licity. However, this occurs only at small radii and does not 
influence average abundance measurements integrated out to 
large radii where most of the cluster mass resides. 

With these cluster elemental abundances, we investigate the 
source of the metals as a mixture of canonical SN la and SN II, 
and propose alternative sources of metals necessary to match 
the observations. 

2. THE ELEMENTS 

The strong n=2 to n=l Ly-a (or K-a) lines of the elements 
from C to Ni lie in the X-ray band between 0.1-lO.OkeV. 
These are the largest equivalent width lines in the X-ray spec- 
trum for clusters with temperatures greater than ~2keV, and 
the most useful for determining elemental abundances. The 
strength of these lines depends on the abundance of the ele- 
ments and their ionization balance, which in turn depends on 
the temperature of the cluster The deep gravitational potential 
well of clusters heats the gas and leaves it highly ionized. The 
gas emits primarily by thermal bremsstrahlung, and for the 
temperature range of galaxy clusters the ionization balance is 
such that most elements have a large population of their atoms 
in the H-like and/or He-like ionization states over most of the 

^ The procee dings of the Ringberg Cluste r Conference, IDeGrandl20O3D 
can be found at: |http : / / www . xray .mpe . mpg . de/ ~ ringberg03/ 



cluster volume. Clusters are optically thin and nearly isother- 
mal, with the result that the line emission is easily interpreted 
without complicating factors such as radiative transport and 
the imprint of non-thermal emission. 

While all the elements from carbon to nickel have their 
main lines in the X-ray band, not all of them are easily visi- 
ble. Elements like fluorine and sodium have abundances more 
than an order of magnitude below the more abundant elements 
such as silicon and sulfur, and are so far not detected in ob- 
servations of galaxy clusters. The list below introduces the 
more abundant and important elements, and the prospects for 
measuring their X-ray lines with ASCA in galaxy clusters. 

2.1. Carbon, Nitrogen and Oxygen 

Low temperature clusters and groups may have nitrogen 
and carbon K-a lines with significant equivalent width. How- 
ever, these lines lie below the usable bandpass of the ASCA 
detectors. 

H-like oxygen has strong lines at 0.65 keV and is an impor- 
tant element for constraining enrichment scenarios because it 
is produced predominantly by type II supernovae. However, 
the response of the ASCA GIS detector is uncertain at these 
energies, and the efficiency of the SIS detector varies with 
time at low energies and is also relatively uncertain. Unfor- 
tunately, the usable bandpass we adopt for ASCA does not go 
low enough to include oxygen. 

2.2. Neon and Magnesium 

The K-ai H-like line for neon is at 1.02 keV and falls right 
in the middle of the iron L-shell complex ranging from about 
0.8-1.4keV. The resolution of ASCA and the close spacing of 
the iron lines makes neon abundance determinations from the 
K-shell unreliable. 

With its K-shell lines also lying in the iron L-shell complex 
(1.47keV), magnesium suffers from the same problems as 
neon and is not well determined with ASCA data. Results for 
both neon and magnesium from the high resolution RGS on 
XMM show that the CCD abundances do not match those ob- 
tained with higher resolution gratings ( Sakelliou et al. 2002^), 
indicating that CCD abundances such as those obtained from 
ASCA are not capable of giving acceptable results. 

2.3. Aluminum 

Aluminum has a higher solar abundance than calcium and 
argon (the two lowest abundance elements considered in this 
paper), but its H-like Kq, line is blended with the much 
stronger silicon He-like Kq line and is not reliably measur- 
able. 

2.4. Silicon and Sulfur 

After iron, the silicon abundance is the next most well- 
determined of all the elements. Its H-like K-ai line at 
2.00 keV and He-like lines at 1.86keV lie in a relatively un- 
crowded part of the spectrum, and silicon's large equivalent 
width leads to a well determined abundance. 

Next to Fe and Si, the high natural abundance of sulfur 
and its position in an uncrowded part of the X-ray spectrum 
make it a well determined element. Its K-ai H-like line is at 
2.62 keV. 

2.5. Argon and Calcium 

The natural abundance of argon is down almost an order 
of magnitude from sulfur, giving it a lower equivalent width. 
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However, the K-ai H-like line at 3.32 keV is in a clear part of 
the spectrum and measurable. 

Calcium is similar to Ar, with a K-ai H-like line at 
4.10keV. 

2.6. Iron and Nickel 

Iron has the strongest set of lines observable in the X-ray 
spectrum. High temperature clusters above 3 keV primarily 
have as their strongest lines the K-a set at about 6.97 and 
6.67 keV for H-like and He-like iron, while lower tempera- 
ture clusters excite the L-shell complex of m any lines between 
about 0.6 and 2.0 keV. Hwan g^etall Sl99^ have shown that 
ASCA determinations of iron abundances from just the L or 
K-shells give consistent results. Iron and nickel are predomi- 
nantly produced by SN la. 

Like iron, nickel also has L-shell lines that lie in the X- 
ray band. But unlike iron, the abundance determinations are 
driven almost entirely by the He-like and H-like K-shell lines 
at 7.77 and S.lOkeV. This is because the abundance of nickel 
is about an order of magnitude less than iron, and the nickel L- 
shell lines are blended with iron's. Nickel abundances using 
the H-like and He-like lines are most reliable for temperatures 
above ~4 keV since there is little excitation of the K-shell line 
below this energy and because the reference data for the L- 
shell lines is not well constrained. 

3. SOLAR ABUNDANCES 

There has been some controversy in the literature as to the 
canonical values to use for the solar elemental abundances. 
The values for the elemental abundances by number that are 
found by spectral fitting to cluster data do not depend on the 
chosen values for the solar abundances. However, for the sake 
of convenience elemental abundances are often reported with 
re spect to the solar v alues. 

'Mushotz kv et al.l J 19961) in their paper report cluster 
abundances with respect to t he photosphe ric values in 
lAnders & Grevesse ( 198^. In lAnders & Gr evesse ( 198^, 
the authors comment on how the photospheric and meteoritic 
values for the solar abundances were coming into agreement 
with better measurement techniques and improved values of 
physical constants, and give numbers for both the photo- 
spheric and meteoritic values. While almost all the elements 
were in good agreement, the iron abundance still showed 
discre pancies between the photos pheric and meteoritic val- 
ues. Ishimaru & Arimoto| (Il997h questioned the claims in 
[Mushotzkv et al. ( 1996) by noticing that they used the photo- 
spheric values when analyzing the data (the default in XSPEC 
then and now), but that the theoretical results they were com- 
paring to used the meteoritic abundances. Since the discrep- 
ancy in the two values for iron was significant, and because 
many of the abundance ratios used in the analysis were with 
respect to iron, the conclusions were based on incompatible 
iron data. 

Since 1989, the situation has improved. Reanalysis of the 
stellar photospheric data for iron that includes lines from Fe II 
in addition to Fe I as well as improved modeling of the solar 
lines ( Grevesse & Sa uval 1999) have bro ught the meteoritic 
and ph otospheric values into agreement. iGrevesse & Sauvall 
llI998h incorporate these changes and others and has become 
the de facto standard for the standard solar composition. Ta- 
bleQ gives the abundances from both sources. 

However, the past history of changes in the adopted solar 
composition implies that there might still be changes in the 
abundance values for some elements. Because of this, and the 



Table 1. Solar Abundances 



hlement 


Anders & Grevesse 
(1989)" 


— 

Grevesse & Sauval 

(1998)'' 


11 


12.00 


12.000 


C 


8.56 


8.520 


N 


8.05 


7.920 





8.93 


8.690 


Ne 


8.09 


8.080 


Mg 


7.58 


7.580 


Si 


7.55 


7.555 


S 


7.21 


7.265 


At 


6.56 


6.400 


Ca 


6.36 


6.355 


Fe 


7.67 


7.500 


Ni 


6.25 


6.250 



Re ferences. — (1) (Anders & Gr evess'3ll989t (2) IGrevesse & Sauvall 
IT99I 

Note. — Abundances are given on a logarithmic scale where H is 12.0. 

"These numbers ai'e the photospheric values, used as the default in XSPEC. 

''These numbers are a straight average of the photospheric and mete- 
oritic values (excep t for oxygen, which has the updated value given in 
lAllende Prieto. Lambert. & AsplunA2001i) . 



problems of comparing results produced with different, in- 
compatible solar values, we quote our results for the elemen- 
tal abundances by number with respect t o hydrogen. We also 
give th e abundances with respect to the lAnders & Grevess3 
Jl989h solar abundances to ease comparisons with previous 
wor ks, and in addition lis t our results with respect to the stan- 
dard Grevesse & Sauval ( 1998) values for convenience and 
for constructing abundance ratios. 

4. OBSERVATIONS AND DATA REDUCTION 
4. 1 . Sample Selection 

We use for our sample all the cluster observations 
in the archives of the ASCA satellite. In iHorner et alJ 
(^pJS submitted)*' (hereafter ACC for ASCA Cluster Cata- 
log), we describe our efforts to prepare a large catalog of 
homogeneously analyzed cluster temperatures, luminosities 
and overall metal abundances from the rev2 processing of 
the ASCA cluster observations. There we give the full details 
of the data selection and reduction; only a brief summary is 
given here. In this paper we use the ACC sample, but our 
focus is the determination of the abundances for individual 
elements in addition to iron. 

The ASCA satellite was launched in February 1993, and 
ceased scientific observations in July 2000. Over the course of 
its lifetime it observed 434 clusters in 564 observations. The 
cluster sample prepared in ACC selects 273 clusters based on 
the suitability of the data for spectral analysis by removing 
clusters with too few photons to form an analyzable spectrum, 
clusters dominated by AGN emission, etc, and is the largest 
catalog of cluster temperatures, luminosities, and abundances. 
However, because the catalog was designed to maximize the 
number of clusters obtained from the ASCA archives, it is 
not necessarily complete to any flux or redshift and could be 
biased because of the particular selections of the individual 
ASCA observers who originally obtained the data. 

The cluster extraction regions in the ACC sample were se- 
lected to contain as much flux as possible in order to best rep- 

* The results in IHomer et alj JAgJS^ubmitte^ are primarily from 
Don Homer's Ph.D. dissertation ^|^fo men2001fl . found online at: 
jhttp ://sol.stsci. edu/ ~horner| 
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Table 2. Stack Parameters 



Stack 


Temperature 


Number of 


Total 


Name 


Bin (keV) 


Clusters 


Counts 


A 


0.5 


17 


50802 


B 


1.5 


44 


228685 


C 


2.5 


35 


261267 


D 


3.5 


47 


478274 


E 


4.5 


38 


277014 


F 


5.5 


37 


391047 


G 


6.5 


39 


261484 


H 


7.5 


20 


111593 


I 


8.5 


14 


93481 


J 


9.5 


13 


171669 


K 


10.5 


22 


135321 



resent the total emission of the cluster. Radial profiles of the 
GIS image were made and the spectral extraction regions ex- 
tended out to the point where the cluster emission was 5a 
times the background level. Standard processing was applied 
to the event files (see ACC). For the GIS detector, we used 
the standard RMF and generated an ARE file for each cluster. 
For the SIS detector, we generated RMFs for each chip and an 
overall ARE for each cluster. Backgrounds for the GIS were 
taken from the HEASARC blank sky fields except for low 
galactic latitude sources (\h\ < 20°) where local backgrounds 
were used. For the SIS, local backgrounds were used unless 
the cluster emission filled the field of view. In ACC, clusters 
with more than one observation had their spectral files com- 
bined before analysis; the joint fitting procedure described be- 
low allows us to deal with multiple observations without a 
problem. 

4.2. Stacking Analysis 

Only the very few brightest cluster observations in the 
ASCA archives have enough signal to noise to allow spectral 
fitting of the intermediate elements. In order to improve our 
sensitivity to these elements, we jointly fit a large number of 
clusters simultaneously. We divide the 353 observations (of 
273 clusters) in the ACC into fitting groups called stacks by 
placing together clusters with similar temperatures and over- 
all abundances (see Figure |2j- Not only does this allow us 
to increase our signal to noise level, but also decreases our 
susceptibility to systematic errors resulting from inaccuracies 
in the instrument calibration. Because clusters with differ- 
ent redshifts (but similar temperatures) are analyzed jointly, 
any energy-localized error in the instrument calibration will 
have less effect on abundance determination because the error 
is not likely to affect all the clusters in a single stack. This 
method of analysis also smoothes over biases that may result 
from the different physical conditions found in clusters {e.g., 
substructure and eccentricity). 

Clusters with more than 40 k counts in GIS2 were not in- 
cluded in our analysis so that very bright clusters do not un- 
duly bias the results for a particular stack. Figure [0 shows 
a histogram of the number of counts per cluster observation, 
and indicates that we only lose 25 data sets by excluding ob- 
servations with more than 40k counts. 

The number of stacks was motivated by our desire to have 
a reasonable number of stacks covering the 1-lOkeV temper- 
ature range, and by the limitations of the XSPEC fitting pro- 
gram (jointly fitting more than about 20 clusters with a vari- 
able abundance model exceeds the number of free parameters 
allowed). We divide the clusters into 22 stacks by first separat- 




2x10^ 4x10"^ 6x10^ Sxio"^ 1x10 



GIS2 Counts 

Fig. 1. — Counts per cluster histogram. We excluded clusters with more 
than 40 k counts in the GIS2 detector from our analysis so that very bright 
clusters do not unduly bias the joint spectral fitting results. The vertical line 
shows the sample cut. There are only 25 cluster observations excluded from 
the analysis because they have too many counts. 



ing them into one keV bins (0-1 keV, 1-2 keV, . . . , 9-lOkeV, 
10+keV), and then dividing each one keV bin into a high and 
low abundance stack using the iron abundances from ACC. 
The split between high and low abundance was made such 
that there are roughly an equal number of clusters in the high 
and low stack for each one keV bin. In the case where there 
are more clusters in a stack than it is possible to jointly fit, 
we divide the stack in two and recombine the results for each 
sub-stack after fitting. Tableland Figure|2]show the number 
of clusters in each stack and the dividing line between high 
and low abundance stacks, as well as how the stacks fill the 
abundance-temperature plane. 

We jointly fit the clusters in each stack to a vari- 
able abundance model modified by galactic absorp- 
tion (tbabs*v apec ( iWilms. Allen. & McCravl 120001: 
ISmith etal]l2001l) within XSPEC). The vapec model uses the 
line lists of the APEC code to generate a plasma model with 
variable abundances for He, C, N, O, Ne, Mg, Al, Si, S, Ar, 
Ca, Fe, and Ni. We fixed He, C, N, O, and Al at their solar 
values and allowed the other elements to vary independently 
(except for Ne and Mg, which we tied together). After 
separately extracting each detector, we combine the two GIS 
data sets and the two SIS data sets together before fitting. 
The redshift for each cluster was fixed to the optical value 
found in the literature, except for the few clusters without 
published optical data (which we allowed to vary [see ACC]). 
The column density was fixed at the galactic value for the 
GIS detectors, but allowed to float for the SIS detectors in 
order to c ompensate for a varying low energy efficiency 
problem (Yaaoobetal. 2000)^. The data for most clusters 
was fit between 0.8-lO.OkeV in the GIS and 0.6-lO.OkeV in 
the SIS. Observations made after 1998 had a higher SIS low 
energy bound of 0.8 keV because of the low energy efficiency 
problem, and about 10 other clusters had modified energy 
ranges because of problems with the particular observation 
(see ACC). 

After fitting each of the 22 stacks, we combined the results 

^ ASCA GOF Calibration Memo (ASCA-CAL-00-06- 
01, vl.O 06/05/0 0) lYaaoob et aU ITOOffl can be found at: 
|http : / /heasarc/ docs/asca/calibration/nhpa ram. html| 
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Fig. 2. — Stack selection diagram. Tlie lines on this abundance-temperature plot show where the boundaries were placed for the individual stacks. Each point 
is a single cluster measurement from ACC. Stacks A through V were fit individually in XSPEC, and the results from the low and high abundance stacks were 
combined into a single result for each one-keV bin (e.g. stack A combined with stack L, etc). 



from the low and high metallicity stacks in the same temper- 
ature range in order to further improve the statistics. The 
difference in the several elemental abundances between the 
low and high metallicity stacks was consistent with an overall 
higher or lower metallicity (ie, the low metallicity stacks L- 
M have slighter lower Fe, Si, S, and Ni than the high metal- 
licity stacks A-K). For example, stack A (the 0-1 kev high 
metallicity stack) was combined with stack L (the 0-1 keV 
low metallicity stack) into a new stack A. These final results 
for 1 1 stacks are given in the next section. 

5. RESULTS FOR INDIVIDUAL ELEMENTS 

We present results for the abundance of the elements Fe, Si, 
S, Ar, Ca, and Ni as a function of cluster temperature. Other 
elements with lines present in cluster X-ray spectra (e.g. Ne, 
Mg and O) have statistical or systematic errors too large to 
allow meaningful results. The main results of our analysis are 
presented in Table |3] which lists the metal abundances of the 
cluster stacks. These results are given by number with respect 
to hydrogen. In Tablel^we give the same results with respect 
to the photospheric solar abundances in lAnders & Grevessd 
{1989.) in order to allow easy comparison with other results 
in the literature. Finally, in Table [S] we give the cluster abun- 



dances with respect to the standard solar comp osition given in 
our Table^adopted from Grevesse & Sauvaltl99 8). 

Figure|3ldisplays the same information as in Table|5l but in 
graphical form. Several interesting points are immediately ap- 
parent in the data. First, the results for iron closely follow pre- 
vious results for cluster metallicities, but have much smaller 
error bars. The actual numerical result for the iron abundance 
is consistent with the previous results of about 1/3 solar, but 
the improved quality of the data allows the detection of trends 
in the iron abundance with temperature: at high temperatures 
above 6 keV, the abundance is constant at a value of 0.3 solar 
Between 2.5 keV and 6 keV, the iron abundance falls from 0.7 
solar to 0.3 solar, and from 0.5 keV to 2.5 keV the abundance 
rises steeply from 0.3 solar to 0.7 solar 

Also of note are the results from the second and third most 
strongly detected elements, silicon and sulfur Here, the abun- 
dance ratios give some important information: the value of 
[Si/Fe] is generally super-solar, while the value of [S/Fe] is 
sub-solar. Table|6lgives the ratios of silicon, sulfur, and nickel 
to iron. Figure |3 shows the trend of the abundance ratios 
[Si/Fe] and [S/Fe] as a function of the iron abundance. Silicon 
and sulfur also show disagreement in their trend with temper- 
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Table 3. Galaxy Cluster Elemental Abundances by Number 



Temperature kT" Silicon'' Sulfur'' Argon'' Calcium'' Iron'' Nickel'' 

Bin (keV) (keV) 



0.5 


n 830-84 
U.S3o82 


95 S'O^-' 


107.9^f,o 


154.8ini 


263 1 305-0 
'^"-'■^224.0 


102 O'O^-l 


l'-3 


1.5 




132 3'38-4 




on 943.2 
-'^■^35.3 


CI q58.0 
5^-'<47.6 


135 6'38-'' 
lJJ.O(324 


2 6"-' 


2.5 




163 9'83.l 


51 6*'''-0 




7 4IO.9 
'•^4.9 


220.822|7 


9.7^3^4 


3.5 


3.683;™ 


206 1^25-3 


42 05' -2 




1 l3-5 


200.7205.8 


20.0?!^ 


4.5 


'*■->' 4.54 


197 6233-1 


34 250-3 
-"^■^20.4 


0.03;2 




127 7'31-4 


I64I9-7 


5.5 


-'•"5.72 


23Q q273.6 
"^•^209.0 


16 232-8 






137.0l«:3 


24.5i:| 


6.5 


7,6.78 


212 5254-O 


26 9'"-'' 
ZD.y23.4 


OO'O 


0'-* 


88.9|2-6 


16.0!!^ 


7.5 


7.45?:g 


1 CO ^229.6 


51 2'*''-2 


U.Uqq 


04-5 


92.ig:i 


23.2?l:i 


8.5 


8.30^:fo 


289.93§8.2 


66.2ii« 


2 025-« 


1 l'0-2 


74 8*2-8 




9.5 


9.63|-80 


364.0^^1^ 




OS''-' 
"■-'0.0 


04-4 


106 2"3-' 


25.8jg7 


10.5 


10 92"-" 


298 n'*05.9 


64.2|i3-0 


0'2-4 


0'-** 

"-'^0.0 




21 328-9 



Note. — The numbers in the sub and superscripts for the abundances are the low and high extent of the 90% confidence region for that element. 
"The values in the temperature column are the fitted temperature of the simultaneously fit clusters in this temperature bin. 
''AH abundances are 1 X lO' times the number of atoms per hydrogen atom. 



Table 4. Classical Galaxy Cluster Elemental Abundances" 



Stack Temperature Silicon" Sulfur" Argon" Calcium" Iron" Nickel" 
Name Bin 



A 


0.5 


270-30 

u.z;(,24 


67O '5 


4 26"* 


11.48»,3i 


220-23 
""1 


010.07 


B 


1.5 


370-39 


440-47 


1.08|!|;„ 


2 302-53 
z.jUj og 


29O 30 
"•^'0.28 


140-23 

"•^^0.02 


C 


2.5 


0.460:52 


320-37 


03'"3 


320-48 


470-49 


540-75 


D 


3.5 


0.580-64 


260-32 


01" 07 


050-'5 

"-"-^0.00 


430-44 


1 13'-32 


E 


4.5 


0.560-$6 


210-31 

U.Zlg ,3 


qOqO 10 


OOQO-'O 
"-""0.00 


270-28 

"•^'0.26 




F 


5.5 


0-6<59 


100-20 


00" 07 


000-08 

"-""0.00 


290-30 

"•^'^0.28 


1.38l;« 


G 


6.5 


600-''2 


170-29 


oo"''' 
"■""0.00 


000-07 
"-""0.00 


I9O 20 


90' 12 
"■^'^'0.70 


H 


7.5 


450-65 


320-52 


00" 23 
"■""0.00 


000-20 
"■""0.00 


2OO 21 


1 31'-'^'* 


I 


8.5 


0.82i;0? 


4,0.71 


oe"^" 

"•"°o.oo 


"■"•'0.00 


160-"* 


1 23 '-*7 
^•"0.79 


J 


9.5 


1 03'-2S 


320-57 


01" 21 
"■"'0.00 


OOQO-" 
"■""0.00 


230-24 


1 45' *5 


K 


10.5 




400-™ 


00" 34 

"■""0.00 


000-25 
"-""0.00 


n 1 80.20 


1 20' ^2 



Note. — The numbers in the sub and superscripts for the abundances are the low and high extent of the 90% confidence region for that element. 
'All abundances are with respect to the solar photosphere elemental abundances given in lAnders & Grevessell989l 



Table 5. Current Galaxy Cluster Elemental Abundances" 



Stack Temperature Silicon" Sulfur" Ai'gon" Calcium" Iron" Nickel" 
Name Bin 



A 


0.5 


270-30 

u.z;q24 


0.59066 


6 1 67 07 

O.lOj 3j 


11 62'3.47 

ll.OZg gg 


320-34 


OlO 07 
"•"^0.00 


B 


1.5 


370-39 




1 561 72 
I.JO, 4, 


9 T-i2.56 
z.jj^ 10 


430.44 


140-23 


C 


2.5 


0.460;51 


0.280;33 


0040"' 

"•"^0.02 


0.320^48 


700.72 


540-75 


D 


3.5 


0.570-63 


0-23^:?^ 


OOlO'O 
"•"^0.00 




63O 65 


1 13'-32 


E 


4.5 


550-'55 


190-27 


000-'5 
"•""0.00 


OQO-'O 
"■""0.00 


400-'*2 


0-92i:^^ 


F 


5.5 


670-76 

"-°'0.58 


ngO is 

"■"'0.07 


OOO-'O 
"■""0.00 


000-08 
"■""0.00 


0.430;« 


i.38|:6^ 


G 


6.5 


590-^' 


0.15026 


000-20 

"•""0.00 


000-07 

"-""0.00 


0.280-29 


90' '2 
"■^"0.70 


H 


7.5 


0.440-64 


0.280-^6 




000-20 
"-""0.00 


U.Zi-Q 27 


1 3l'-68 


I 


8.5 


0<fe 


0.36062 




050-^5 

"-"•^0.00 


240.26 


1 23'-67 
^•^•'0.79 


J 


9.5 


l-01a?9 


0.280-50 


010 30 
"•"^0.00 


000-'9 
"-""0.00 


340-36 


1 45'-'*5 

^•^-^L05 


K 


10.5 


0.83ij3 


0.350:61 


OOO 49 
"■""0.00 


000-25 
"-""0.00 


270-30 

U.Z/Q24 


1 20' -62 



Note. — The numbers in the sub and superscripts for the abundances are the low and high extent of the 90% confidence region for that element. 
"All abundances are with respect to the average of the photospheric and meteoritic solar elemental abundances given in TablelTladapted from iGrevesse & Sauval 
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Table 6. Elemental Abundance Ratios 



Stack Temperature [Si/Fe] [S/Fe] [Si/S] [Ni/Fe] [Fe/H] 

Name Bin 



A 


0.5 


-0.08:003 


26"-32 


_0 34-0.27 


-1.81^^64 


-0 49-" 47 


B 


1.5 


-0.07^:°^ 


A A4-O.OI 


-0 02" "' 

"■''^-0.06 


_A 47-0.26 


-0 17-".36 


C 


2.5 


-0.18^:^3 


-0 40-" 33 


u.zig 12 


-0 11 ""3 
"■''-0.23 


-0 le-"'* 

"■'"-0.17 


D 


3.5 


-0 04-" '"' 


-n 44-"-36 
"■^^-0.55 




25"-32 


-0 20-"" 
"■^"-0.21 


E 


4.5 


13"-' 


-0 M-°-" 


47"-65 




-0 39-"-38 


F 


5.5 


19"-25 
"■'^0.13 


-o.69:«:3^ 




50" " 


"■^ -0.38 


G 


6.5 


32"-''" 




O61O.86 

U.OIqji 


50"-^" 


-0 55-"-^3 


H 


7.5 


A 10O.34 
O.IS_(,o2 


-0.02"„2§, 


20O 45 

-0.29 


65""' 


A 54-0.50 

"■^^-0.57 


I 


8.5 




0.18"„«3 


35"''' 


71" 


A 63-0.58 
"■°-'-0.68 


J 


9.5 




-0 07" '* 


55"-*' 




A 47-0.44 

"■^'-0.51 


K 


10.5 




1 l"-36 
u. 11 


0.38"„«t2 




-0 57-"-53 
"■^'-0.62 



Note. — All abundance ratios ai'e with respect to the cuirent abundances given in Table Isl The numbers in the sub and superscripts for the abundances 
are the low and high extent of the 90% confidence region for that element. Abundances are given in the usual dex notation, ie: [A/B]= logio(NA/NB)ciustcr ~ 
logio(NA/NB)0. 
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Fig. 3. — The galax y cluster elemental abundances as a function of temperature. These abundances are with respect to the solar abundances of 
IGrevesse & Sauval 119981) . The eiTor bars are the 90% confidence interval for that elemental abundance; the eiTor bars for iron aie smaller than the plotted 
points. 



ature — tlie relative abundance of silicon rises with tempera- 
ture, and sulfur falls. We find that the silicon abundance rises 
from '^0.3 solar in cooler clusters to ^0.7 s olar in hot clus- 
ters, in excellent agreement with the results of lFukazawa et alJ 



(1998). 

Calcium and argon are noticeable for their lack of a de- 
tectable signal. Both are not detected in the stacked ASCA 
data, yielding only upper limits over the temperature range 
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Fig. 4. — The silicon and sulfur abundances ratios with respect to iron. 



from 2-12 keV. 

The relative abundances for nickel are measured to be 
higher than any other element in our study. For clusters 
hotter than 4keV, nickel is about 1.2 times the solar value. 
Lower temperature clusters are too cool to significantly ex- 
cite the nickel K-a transition, and abundances derived from 
the Ni L-shell lines are not reliable. The Ni/Fe values are ex- 
tremely high at about 3 times the solar value, confirming the 
results of Dupke & White (2000a) : Dupke & Arnaud (2001) : 
lT5^nke& White! Il2000bft. 

6. SYSTEMATIC ERRORS 

The abundances results in this paper depend upon the mea- 
surement of particular spectral lines that are not always very 
strong. Because of this, a proper understanding of the results 
must take into account any systematic errors that may bias 
the results. Of particular concern to this work are systematic 
errors in the calibration of the effective area, since these of- 
ten manifest themselves as lines in residual spectra of sources 
with a smooth continuum. If these residual lines fall at the 
same energy as the important elemental spectral lines they can 
have a significant affect on the derived elemental abundances. 

The ACC catalog paper discusses general ASCA systemat- 
ics and our sensitivity to them. In addition, there are several 
other tests we have undertaken to determine the effect of small 
line-like systematic errors in the effective area. In order to de- 
termine the size of any calibration errors, we have fit spectra 
from broad band continuum sources and quantified the resid- 
uals. We have also used these residuals to correct the cluster 
spectra, and have compared the derived elemental abundances 
to those from the uncorrected spectra. 

We have used an ASCA observation of Cygnus X- 1 (a bright 
source where the systematic errors dominate the statistical 
ones) to measure the size of residual lines in fitted spec- 
tra of a continuum source. These lines could be due to er- 
rors in the instrument calibration and could affect the abun- 
dance determinations in clusters. Using a power law + 
diskline model in XPSEC, we measure the largest residual 
in the Cyg X-1 spectrum to have an equivalent width of 17 eV 
at 3.6 eV. Using the Raymond-Smith plasma code to model 
the equivalent width of elemental X-ray lines, we find that the 
only element with a small enough equivalent width to be pos- 



FlG. 5. — An absorbed power law fit to the ASCA 3C273 data. The data 
and model are .shown in the top panel for the 4 ASCA detectors, and the ratio 
of the data to the model is shown in the bottom panel. The lack of any signif- 
icant residuals indicate that the ASCA effective area calibration is free of any 
significant line-like systematic errors. 



sibly affected by a 17 eV residual is calcium, with an equiva- 
lent width of 25 eV for very hot clusters (> 6keV). However, 
the 17 eV residual lies at the wrong energy to affect calcium 
(lines at 3.8, 3.9, and 4.1 keV). In addition, the positive resid- 
ual would serve to increase the measured calcium abundances; 
our measured calcium abundances are lower than expected. A 
similar test with the bright continuum source Mkn 421 (with 
the core emission removed to prevent pileup) shows a maxi- 
mum residual of 18 eV equivalent width at the 2.1 1 keV Au 
edge of the mirror These results indicate that any line-like 
calibration errors are not large enough to affect the cluster el- 
emental abundance measurements in this study. 

We have also used an ASCA observation of 3C273 (se- 
quence number 12601000) as a broad band continuum source 
to check for calibration errors in the response matrices and 
their effect on the derived cluster elemental abundances. We 
fit 3C273 with an absorbed power law model (Figure|5} and 
extracted the ratio residuals to the best fit. We then used these 
residuals as a correction to the cluster data, and then refit the 
clusters to find the abundances. The abundances from the 
corrected data are completely consistent with the abundances 
from the uncorrected data, further indicating that errors in the 
effective area calibration do not affect the derived elemental 
abundances. 

We have also investigated the abundances derived using the 
GIS and SIS detectors separately in order to check the consis- 
tency of our results. While the SIS has better spectral resolu- 
tion, it also suffers from slight CTl and low energy absorption 
problems that do not affect the GIS. Figure|6lshows these re- 
sults. In general, the GIS and SIS abundances match very well 
for most of the elements. However, nickel and silicon show a 
systematic trend of slightly higher SIS abundances. The SIS 
nickel abundances are not as reliable as the GIS abundances 
because the GIS has more effective area than the SIS at Ni K. 

The systematic trend in the medium te mperature sihcon 
abundances is more difficult to understand. 'Fukazawa' ( Il997h 
showed that the GIS and SIS sulfur and silicon abundances 
for his cluster sample were well matched. Individual analy- 
sis of the bright, medium temperature clusters Abell 496 and 
Abell 2199 show that the results are indeed real and verify that 
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Fig. 6. — A comparison of abundance results using the GIS alone, the SIS alone, and the GIS and SIS combined. The GIS and SIS agree well except for nickel 
and a slight trend in the medium temperature silicon results. The nickel discrepancies are due to the lower effectiveness of the SIS at Ni K, and the SIS silicon 
measurements suffer from proximity to the detector Si edge. The errorbars for the GIS and SIS alone points are omitted for clarity, but are about a factor of \/2 
larger than those for the combined points. 



the SIS gives higher Si abundances than the GIS. Additional tion to this problem is that the medium temperature clusters 
fits to the data that allow the SIS gain to vary also do not sig- most affected by this trend mostly lie at a redshift such that the 
nificantly affect the abundance results. One possible contribu- Si K-shell lines (1.86 and 2.00 keV) are redshifted very close 
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to the Si edge (1.84keV) in the detector. These silicon abun- 
dances are less reliable because of structure prese nt in the Si 
edge not well modeled in the instrument response jMori et alJ 
H)00). 

The fit results in the 0.5 and 1.5keV bins from Tables|3]S 
and |5] show elevated calcium and argon abundances. These 
very high fit results (1-10 times solar) are not believed to 
be indicative of the actual cluster abundances because of sys- 
tematic errors that preferentially affect these low temperature 
bins. Two sources of systematic error we have identified both 
become important at low cluster temperatures and affect the 
modeled spectrum at the higher energies of the Ar and Ca K- 
shell lines. 

The first is background subtraction: because the dominant 
bremsstrahlung emission is limited to lower energies for the 
cooler clusters, the cluster flux present at the higher energies 
of the Ar and Ca K-shell lines becomes comparable to the 
background. Small errors in the blank sky background can 
then significantly affect the background-subtracted data, lead- 
ing to incorrect abundances at higher energies where the rela- 
tive errors from the background subtraction become large. 

Also important at the lower cluster temperatures 
where the bremsstrahlung emission does not domi- 
nate throughout the X-ray spectrum is the contribu- 
tion from point sources within the individual galax- 
ies. [ Angel ini. Loewenstein. & Mushotzkv (2001) and 
llrwin. Athev. & Bregman li20()3l) have shown with data 
from elliptical galaxies that the contribution from X-ray 
binaries is important and can be characterized with a 
7keV bremsstrahlung model. Since we do not include this 
component in our cluster fits, we expect that our Ar and 
Ca abundance for low temperature clusters will be driven 
upwards in an attempt to try and match the flux actually 
contributed by the X-ray point sources. 

7. COMPARISONS 

We believe that X-ray determinations of elemental abun- 
dances in galaxy clusters are one of the most accessible means 
for obtaining useful abundance measurements. While our re- 
sults are consistent with and significantly expand previous 
cluster X-ray measurements, they do not always agree with 
measurements of elemental abundances in other objects. 

7.1. Comparison to Other Cluster X-ray Measurements 

7.1.1. ASCA Measurements 

Measurements of abundances for elements other than iron 
has historically been difficult because of the low equivalent 
width of the lines and the limited resolution of X-ray detec- 
tors. Aside from a few limited results with crystal spectrom- 
eters on Einstein and other satellites, the first chemical abun- 
dance measurements of el ements oth er than iron in galaxy 
clusters was made bv (Mushotzkv et al.l ( 1996) and made use 
of the high resolution of the CCD cameras onboard ASCA. 
Their results for 4 moderate temperature (2.5-5. OkeV), very 
bright clusters show an average Fe abundance of 0.32 solar, 
0.65 solar for Si, 0.25 for S, and 1.0 for Ni. AH of these re- 
sults are in agreement with our results for moderate tempera- 
ture clusters presented in Table|4] F ukazawa et al. (.199 8) also 
reported results for the silicon abundance in clusters that used 
ASCA data. Their results for 40 clusters was in general agree- 
ment with the data from Mushotzkv et al. ( 1996), showing a 
generally constant iron abundance for clusters above 3keV, 
and also hinted at the silicon trend with temperature presented 
with more detail here. 



IFmosuenov, David. & PonmanI ( l2000h and 

iFinoguenov. Arnaud. & D avidI (120011) concentrated on 
the detection of abundance gradients with ASCA. Their data 
for Fe and Si at large radii are in agreement with ours. 
However, they present measurements of S coupled with Ar 
in order to reduce errors in the S determination. We see that 
these elements have very different abundances, and can't 
compare our data directly to theirs. 

7.1.2. XMM Measurements 

Further results for the intermediate element abundances in 
clusters with ASCA were hampered by its moderate resolu- 
tion, and especially because only a handful of clusters have 
high enough X-ray flux to allow meaningful measurements. 
The advent of Chandra and especially XMM have changed 
this with their higher resolution CCD cameras and much in- 
creased effective area. While measurements of argon and cal- 
cium are still beset with systematic problems in the XMM re- 
sponse and background subtraction, measurements of silicon 
and sulfur are largely free of these complications. 

The much improved spatial and spectral resolution of XMM 
allow not only abundance measurements averaged over the 
whole cluster, but also spatially resolved ones. While our 
measurements in clusters are overall spatial averages, the 
measurements with XMM are usually spatially resolved and 
give information on abundance gradients within clusters. 
While this information is important and sheds valuable light 
on the enrichment mechanisms in clusters, we will limit our 
comparisons to overall abundances measured with XMM, or 
with abundances in the more voluminous outer regions of 
clusters with spatially resolved abundances. 

The many XMM pap ers on M87 ( Bohringer et al. 200 ll 
Molend i & Gastaldello"2001': 'Gastaldello & Molendi, .20o3 
Matsus hita. Finoguenov. & Bohringer 2003) all find a signif- 
icant abundance gradient in the center of the cluster In 
the very center, Molendi & Gastaldello (2001) find generally 
higher abundances than in the outer regions, with the iron 
abundance at about 0.6 the solar value. They find silicon to be 
1.0 in the same region, giving [Si/Fe] only slightly lower than 
our results for a similar temperature cluster. However, sulfur 
is about 1 . 1 solar, overly abundant with respect to iron than in 
our results. The discrepancies between some of our results for 
2.5 keV clusters and the XMM M87 results are tempered by 
the fact that earher ASCA results for M87 (jMatsumoto et al] 
'1996VGuainazzi & Molendi'1999^ are similar to the XMM re- 
sults, and by the fact that the M87 data is for only one cluster, 
while our ASCA results are for an average of several clusters 
at the s ame ternperatur e. 

Tam ura et alJ ( 120011) find sulfur and silicon to have similar 
abundances in the center of A496 with RGS data from XMM, 
slightly different than our results that show silicon slightly 
higher than sulfur by a factor of about 1.5. Unfortunately, the 
M87 data concentrates on the cluster core, and the A496 data 
is taken with the RGS, which also is limited to the centers 
of bright clusters. This bias towards the very centers of clus- 
ters makes a comparison with our average cluster abundances 
difficult 

Fino guenov et alJ (120021) also look at XMM data for M87, 
but focus their paper on a discussion of the supernovae enrich- 
ment needed to account for the observed abundances. Their 
data leads them to conclude, as we do, that the standard yields 
of the canonical SN la and SN II models are not sufficient to 
explain the pattern of abundances observed. Their solution 
calls for an additional source of metals from a new class of 
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Fig. 7. — The silicon and sulfur abundances as a function of temperature, 
compared to their abundance in stars. The upper grey bar is the S/Fe data 
fromlTimmes. Wooslev. & Weaver ( 1995), and the overlapping lower grey 
bar is the Si/Fe data from the same source. 



type I supernovae. This scenario does not fit our data well; 
specifically, it leads to an overproduction of sulfur, argon and 
calcium, with not enough silicon produced to match our ob- 
servations. Finoguenov et al. (2002) focus on the inner the 
inner 70kpc region of M87 because of the cluster's proxim- 
ity. Their results in this region are where the influence of 
the galaxy's stellar population and associated SNla are most 
keenly apparent. Even at the outer parts of this region, the 
abundances show signs of a transition to a more typical clus- 
ter pattern. This inner region contributes a small fraction of 
the emission measure of most systems in our ACC sample, 
and makes comparisons difficult since we sample different ar- 
eas of the cluster 

7.2. Comparison to Measurements in Dijferent Types of 
Objects 

7.2.1. The Thin Disk 

The standard solar elemental composition is based on mea- 
surements taken of the sun, which resides in the thin disk of 
the galaxy. Our cluster results show significant differences 
with stellar data from the thin disk. Figure0shows silicon and 
sulfur abundance ratios in clusters c ompared with data from 
iTimmes. Wooslev. & WeaveJ il995l) . The stellar data is over- 
plotted on the cluster points as grey bars. The [Si/Fe] clus- 
ter data overlaps with the stellar data, but has a much greater 
range. The [S/Fe] cluster data lies almost totally below the 
stellar data. The largest discrepancy is between the cluster 
and stellar nickel abundance ratios; clusters are higher than 
the stellar data of 0.1 by 0.5 dex. This very high value for 
[Ni/Fe] is not seen in stars at any metallicity. 

7.2.2. The Thick Disk 

The majority of the stellar mass in clusters resides in 
elliptical galax ies and the bulge s and halos of spirals 
(iBell, Mcintosh . Katz. & W einberg 12003ft . Prochas ka et alJ 
( l2000i) find that the abundance patterns of the thick disk are in 
excellent agreement with those in the bulge, suggesting that 
the two formed from the same reservoir of gas. If this is 
true, and if the majority of the enriched gas in the ICM origi- 
nated in these numerically dominant stellar populations, then 



we would expect that the abundance pattern in the ICM is in 
good agreement with that in the thick disk of the Milky Way. 
However, reality is more complicated. While Pro chaska et alJ 
(|2000) find that the thick disk has enhanced abundances for 
the a. elements compared to the thin disk and solar neigh- 
borhood, they also find [Ni/Fe] in the thick disk is similar to 
its solar value. Additionally, the calcium abundance ratio is 
found to be super solar ([Ca/Fe] = 0.2), and sulfur is more en- 
hanced with respect to iro n in the thick disk than in c l usters . 
Also, the recent work of Pompeia, Barb uv. & Gren on ( 2003|) 
shows a [Si/Fe] value of nearly 0.0, much lower than what we 
observe here. 

7.2.3. H U Regions and Planetary Nebulae 

Comparison of our cluster results with H II regions or plan- 
etary nebulae is difficult because of the strong effects of dust 
and non-LTE conditions on ab undance measurements of el- 
ements like iron and silicon (Peimbert, Carigi, & Peimberj 
i2001; Stasinska 2002). C, N, and O are more readily observed 
in these objects, but are not observed with ASCA. 

7.2.4. Damped Ly-a Absorbers 

'Prochaska & Wolfe ( 2002) show abundance measurements 
from a large database of damped Ly-a (DLA) observations. 
Their results are for the redshift range 1 .5 < z < 3.5 and indi- 
cate that there is already significant a element enhancement at 
high redshifts. This conclusion is consistent with that reached 
by Mushotzkv et al. ( 1996) in an anal ysis of high redsh ift 
clusters. The [Fe/H] value observed by Prochaska & Wolf^ of 
-1.6 is constant over a wide range of redshift, indicating that 
the enrichment was at earlier times. This value is much lower 
than our value of [Fe/H] = -0.5 (for our higher temperature 
bins where the iron abundance is fairly constant), showing 
the importance of supernovae enrichment in the ICM. Their 
average [Si/Fe] value of 0.35 agrees with our results for mod- 
erate temperature clusters, and their spread of [Si/Fe] mea- 
surements from 0.0 to 0.5 is well matched with ours. Their 
[Si/Fe] values are also constant over a wide redshift range, 
further supporting early enrichment. However, their results 
for sulfur differ from ours in that they show more enrichment 
with respect to iron ([S/Fe] ^ 0.4) than we do. This sug- 
gests that later enrichment by supernovae into the ICM had a 
reduced role for sulfur in comparison to the other elements. 
Nickel is also different in these systems than in clusters, with 
[Ni/Fe] of only 0.07; much different than o ur very high value 
of -0.5. While'Proc haska & Wolfd ( 1200? ) have some results 
for argon, they are widely scattered and not easily interpreted. 
The differences with the cluster measurements of S and Ni in- 
dicate that the stellar population that enriched these systems 
is probably not the origin of the metals in the ICM. 

7.2.5. Lyman Break Galaxies 

iPettini et alJ ( l2002h present detailed observations of a 
lensed Lyman break galaxy at z = 2.7. Their results also in- 
dicate that significant metal enrichment have already taken 
place at high redshift. Their results for O, Mg, Si, S, and P 
are all at about 0.4 the solar value, showing the results of fast 
supernovae processing. However, the iron peak elements are 
not as enriched, with the values of Mn, Fe, and Ni at only 1/3 
solar. The interpretation is that this galaxy has been caught in 
the middle of turning its gas into stars, and that while enrich- 
ment from SN 11 has occurred, not enough time has passed 
to allow significant enrichment from the iron peak producing 
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Fig. 8. — The Type la and Type II supernovae yields. The horizontal axis 
is atomic number — O is 8, Fe 26, Si 14, S 16, Ni 28. The yields are in solar 
masses per supernova. Iron comes mostly from SN la and nickel even more 
so. Si, S, Ar, and Ca are a mix of SNII and SNIa but with a majority SNII 
contribution. The SN la data comes from the W7 model, and the SN II data 
is from the IMF mass-averaged TNH-40 model. 



Fig. 9. — The supemovae-type ratio derived from the silicon and s ulfur 
abundances. The SNIa yields of the W7 model fNomoto et al."1997b^ and 
the SN II yields of the IMF mass-averag ed TNH-40 iNomoto et al. 199?!) 
model were used to compute the Si/Fe and S/Fe abundance ratios for 100 
mixtures ranging from pure SN la enrichment to pure SN II enrichment. The 
measured data was compared to the model output in order to determine the 
relative ratio of SN la and SN 11. A value of 1 on the plot indicates enrichment 
by solely SN la while a value of indicates enrichment solely by SN II. 



SNIa. The overall metallicity of [Fe/H] = -1.2 is closer to 
our observed cluster value than the measurements from the 
DLAs are. However, the abundance ratios in this object are 
also sufficiently different than in clusters to indicate that the 
stellar population that enriched the Lyman break galaxies is 
also probably not the source of the metals in the ICM. 

7.2.6. The Ly-a Forest 

While it is possible to probe the metallicity history 
of the universe as a fu nction of redshift (Songaila 2001; 
ISongaila & Cowijl2002h . the difficulty inherent in measur- 
ing lines from many different elements in the Ly-a forest has 
made abundance measurements of most elements besides Fe, 
C, N, and Si unavailable. Ly-a forest measurements show 
metallicities of [Fe/H] ranging from -2.65 at a redshift of 5.3 
to -1.25 at redshift 3.8. All of these are almost two orders of 
magnitude lower than our measurements. 

8. SUPERNOVAE TYPE DECOMPOSITION 

With the results of §|5]for the abundance of the intermedi- 
ate elements in clusters, it is possible to try to constrain the 
mix of supernovae types that have enriched the ICM. In § 18.11 
we check to see how well the yields from the standard SN la 
and SNII models can reproduce the cluster observations. We 
find that other sources of metals are necessary. Then we com- 
ment on the expected cluster abundances derived using some 
of the standard models in the literature. In § 18.31 we investi- 
gate different SN models that produce intermediate elements 
with the abundance ratios necessary for reconciling the mod- 
els with the observations. 

8.1. SN Fraction Analysis using Canonical SNIa and SNII 

Models 

We use the yields of the revised W7 model ('Nomoto et al.' 
[T^97b) for SNIa and the TNH-40 yields (Nomotoetal. 
[ip97a) (mass-averaged by integrating across the IMF to an 
upper limit of 40 solar masses) for SN 11 (shown in Figure |8j 
as our basis. Each model includes the amount of each of the 



intermediate elements produced in a supernova of that type. 

We used the data from the models to produce a table with 
the yields and abundance ratios for 100 different mixtures of 
SN la and SN 11 ranging from pure SN la enrichment to pure 
SN 11 enrichment. We then compare our two best-measured 
cluster abundance ratios (Si/Fe and S/Fe) with the model ra- 
tios in order to arrive at a SN type fraction. The results are 
plotted in Figure|9]and show the SN type fraction derived from 
both Si/Fe and S/Fe as a function of cluster temperature. 

The difficulty of calculating SN yields leads to acknowl- 
edged unce rtainties for some of the elements of ab out a fac- 
tor of two (Gibs on. Loewenstein. & Mushotzkvl 19971) . In or- 
der to investigate the effects of different SN models on the 
derived SN fraction, we have compared the SN fractions de- 
rived from several different SN models in Figure^] Each SN 
fraction determination uses the W7 SN la model as revised in 
INomoto et aP (Il997bl) . but uses a different SN II model. The 
result of this comparison indicates that different SN II models 
can change the average value for the derived SN fraction, but 
does not change the fact that none of the seven SN II models 
results in a consistent SN fraction from both the cluster Si/Fe 
and S/Fe data. 

The results we have found for the cluster intermediate ele- 
ment abundances do not fit the standard model where most of 
the elements are produced in a simple mix of standard SN la 
and SNII events. The observed silicon abundance is too high 
with respect to iron, and the observed sulfur abundance too 
low. Calcium and argon are also too low. The SN type frac- 
tions derived from the two abundance ratios are not consis- 
tent with each other, and change dramatically as a function of 
cluster temperature. These results indicate that a unique, con- 
sistent decomposition of the ICM enrichment into SNIa and 
SNII contributions is not possible. 

Faced with these inconsistencies, new models must be ex- 
amined that can better explain our results. Some possibilities 
are included in the following section. 
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Fig. 10. — The supemovae-type ratio derived from the silicon and sulfur abundances, comparing the results from many different SNII models. This figure is 
similar to Figurel9l but shows that many different SN II models all have difficulty producing a consistent type fraction from the silicon and sulfur data. Also, the 
trend of cooler clusters predominantly enriched by type SN la products and hotter clusters enriched by SN 11 products does not depend on the SN II model used. 
Again, a value of 1 on the plot indicates enrichment by solely SN la while a value of indicates enrichment solely by SN 11. 
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8.2. Expected Abundances from Standard SN Models 

There are many different supernovae yield computations in 
the scientific literature. For SNll, there are variations in the 
progenitor masses, metallicities and internal structure; in ex- 
plosion energy and placement of the mass cut; and in adopted 
nuclear reaction rates. For SN la, parameters include the cen- 
tral white dwarf progenitor density and burning front propa- 
gation speed. Initially, we consider the standard W7 defla- 
gration SN la model and_Salpeter- IMF-averaged SN 11 yields 
from|Nomoto e t all Jl997bh . In all cases, a combination of 
SN la and SN 11 yields is necessary to match the observed re- 
sults. 

About 0.02 SN 11 events per current solar mass of stars will 
occur if most of the stars in cluster galaxies form over an in- 
terval much shorter than the age of the universe (indicated 
by the predominance of early-type galaxies in clusters) and if 
star formation proceeds with a standard IMF (Kroupa 2002). 
This calculation takes the mass lost to winds and remnants 
into account and assumes the SN 11 progenitors are stars with 
original mass > 8 Mq. Column 2 of Table0shows the obser- 
vational range we aim to explain. Column 3 of Table0shows 
the model abundances from ICM enrichment by SN 11 (for the 
gas-to-star mass ratio of 10 that is typical of a rich cluster). 

If the assumption is made that SN la provide the remainder 
of the iron (The enrichment levels due solely to SN la enrich- 
ment are in column 4 of TableQl, then the total expected abun- 
dances (SN 11 + SN la) are as shown in column 5 of Table 
The number of SN la required corresponds to an average of 
1.8 SNU over 10'° years for a typical ratio of ICM mass to 
stellar blue light of 4 0, significantly higher th an the estimated 
optical rate at z = dCappellaro et al.lll997^ . This combina- 
tion of SN 11 and SN la accounts for the observed abundances 
of nickel and sulfur (compare columns 2 and 5 of Table 0. 
However, only about half of the observed silicon is accounted 
for and calcium and argon are slightly overproduced. Differ- 
ent models must be investigated to make up this significant 
shortfall. The most variation in SN yields is in SN 11 models, 
and we explore some of these alternatives below. 

Unfortunately, many SN 11 models, in par- 

ticular those of [W ooslev & Weaver and 
iRauscher. Heger. Hoffman. & Wooslev. (.20021) . generally 
derive higher yields of sulfur, argon and calcium that ex- 
acerbate the conflict between theory and observations. In 
addition, SN 11 iron yields are uncertain by at least a factor 
of two ( Gibso n. Loewenstein. & Mushotzkv 1997) because 
of their sensitivity to the assumed mass cut. If one arbitrarily 
and exclusively increases the SNll contribution to the iron 
yield by a factor of two, then the SN la contribution to sulfur, 
argon and calcium goes down by a factor of two and is in 
marginal agreement with the observations. Unfortunately, 
this scenario increases the silicon deficit and the nickel 
abundance falls to the unacceptably low value of about half 
solar. 

Models that expUcitly s ynthesize more iron {e.g., those of 
IWooslev & Weav o'^ ( 1995) with enhanced explosion energy) 
also generally produce more sulfur, argon and calcium that 
enhances the probl e m wit h those elements. However, the 
IWooslev & WeaveJ J 19951) enhanced energy models having 
zero metalUcity progenitors have a very different nucleosyn- 
thetic profile — the production of iron and nickel averaged 
over the IMF is doubled without an increase in sulfur, argon 
and calcium. Unfortunately, silicon and nickel remain under- 
produced by a factor of two if these SN 11 yields are adopted. 



Other variations on the SN 11 yields also fail to explain the 
low ratio of sulfur, argon and calcium to silicon. Neither 
adopting delayed detonation SN la models JFinoguenov et all 
2002) nor in creasing the silicon abundan ce by assuming a flat 
IMF (Loewenstein & Mushotzkv 1996) solves the problem. 
Finoguenov et al. ( 2002) explain the abundance pattern and its 
radial variation in M87 by proposing (1) a radially increasing 
SNll/SNla ratio, (2) high Si and S yields from SNla (favoring 
delayed detonation models) and an ad hoc reduction in SNll S 
yields, and (3) a radial variation in SNla yields (corresponding 
to delayed detonation models with different deflagration-to- 
detonation transition densities). While this scenario appears 
promising with regard to the galaxy-dominated inner regions 
of rich clusters an d, perhaps groups, there are difficulties - as 
iFinoguenov et al .' ( 2002 ) acknowledge - in explaining the pat- 
tern in the large scale abundances of high-temperature clus- 
ters. 

One possible solution is to assume a flat IMF to boost SN 11 
silicon production (and account for at least half of the iron), 
combined with ad hoc reductions in sulfur, argon and cal- 
cium yields by a factor of two. An increase in the SNll 
nickel yields (uncertain by a factor of t wo because of mass cut 
and core neutron excess uncertainties. iNakamura et al.ll999h 
would also have to be instated to account for the lowered 
nickel from the reduced SN la output. On the other hand, if a 
standard IMF is maintained, then less significant decreases in 
SN 11 production of sulfur, argon and calcium are needed. For 
this case, the silicon deficit would have to be made up by an- 
other mechanism; enhanced hydrostatic production of silicon 
or a suppressed efficiency of silicon burning are possibilities. 

8.3. Alternate Enrichment Scenarios 

If the observed cluster iron abundance and stellar elemental 
abundances are used to set the contribution from supernovae, 
then the theoretical yields will meet or slightly exceed the 
ASCA abundance observations for argon and calcium while 
silicon itself will be significantly underproduced by a factor 
of ^2. If we accept the standard SNll enrichment given in 
column 3 of Table then we must have an additional nucle- 
osynthetic source that produces substantial silicon, but little 
or no sulfur, argon, or calcium. Some possibilities exist in 
the literature, and we explore some of these scenarios below. 
We have chosen single mass SN II models that can selectively 
enhance silicon, and collectively call these models "SN IIx" 
to distinguish then from the canonical SN II models. For the 
models presented below, we calculate the number of SN IIx 
events necessary to produce the observed amounts of iron and 
silicon that are not produced by regular SN II. 

8.3.1. Very Massive Stars 

iThielemann. Nomoto. & Hashimo to ( 1996) have calculated 
the yields for a 70 Mq progenitor with solar metallicity. The 
abundance ratios in the ejecta with respect to silicon are: 0.4, 
0.28, 0.26, 0.048, and 0.12 relative to solar for S, Ar, Ca, Fe, 
and Ni. With these results, the observational overabundance 
of silicon can be accounted for with ^ 2 x 10""* SN IIx per so- 
lar mass in the ICM. This rate will provide the extra silicon 
necessary, is only about 10% of the total number of SN II ex- 
pected by integrating a normal IMF, and only produces minor 
perturbations in the other elements. The abundances expected 
from the combination of SN la, canonical SN II, and SN IIx 
from very massive progenitors is given in column 6 of Ta- 
bleEl 



1 6 BAUMGARTNER ET AL. 



Table 7. Observed and Model Abundances 



Element 


ASCA 


SNII 


SNIa 


SNII+SNIa 


+SNIIx" 


+SNIIx'' 


+SNIIx'^ 


Si 


0.6^0.8 


0.34 


0.095 


0.4 


0.7 


0.7 


0.7 


S 


0.1-0.4 


0.20 


0.09 


0.3 


0.4 


0.3 


0.4 


Ar 


0.0-0.2 


0.17 


0.075 


0.25 


0.3 


0.2 


0.3 


Ca 


0.0-0.2 


0.19 


0.09 


0.3 


0.35 


0.3 


0.35 


Fe 


0.4 


0.14 


0.26 


0.4 


0.4 


0.4 


0.4 


Ni 


0.8-1.5 


0.16 


0.85 


1.0 


1.0 


1.0 


1.0 



Note. — Abundances are given with respect to the solar values in TableQ column 3. 
"70 Mq progenitors IThielemann. N omoto. & Hashimoto 1996) in addition to SNII and SNIa. 
''0.01 solar abundance, > 30 Mq progenitors iWooslev & Weaver 1995) in addition to SNII and SNIa. 
'^70 M0 He core Pop III progenitors IHeger & Woosle\l2002) in addition to SNII and SNIa. 



8.3.2. Massive, Metal-poor Stars 

Very small amounts of the elements heavier than silicon are 
produced in the high mass (> 25 M(7)), l ow-metallicity (0.01, 
0.1 solar) models of Wooslev & Weaver (1995). These yields 
give less than 10% of the solar ratios (with respect to silicon) 
for elements heavier than silicon. Hence, silicon abundances 
can be increased without causing the other elements to exceed 
their observed values. In order to produce the observed abun- 
dances, about 5.7 X 10"'* SN llx per solar mass in the ICM are 
necessary, assuming progenitors of 30-40 M0 with 0.01 so- 
lar abundance and an IMF that goes as M"' . The abundances 
expected from the combination of SN la, canonical SN II, and 
SN IIx from massive, metal-poor progenitors is given in col- 
umn 7 of Table 

8.3.3. Population in Stars 

The ejecta of supernovae preceded by very massive, zero 
metallicity progenitor stars (Pop III stars) have a substan- 
tially different composition than standard SNII ejecta (e.g., 
Eeger & Wooslev (2002)). If the Pop III mass function is 
dominated by the low m ass (He cores < 80 M©) models of 
iHeger & Wo^^ UM^ . then the relative Si overabundance 
from observations can be explained. If higher core masses are 
used, the silicon overabundances are not as well explained, 
but more modest underabundances are found for sulfur, ar- 
gon, and calcium. For example, about 2.3 x 10"^ SNIIx 
events per solar mass in the ICM from progenitors with 70 
Mq He cores matches our observed abundances and produce 
the abundances given in column 8 of Table[7] 

Figure[^is based on TableQand shows graphically what 
proportion of each element's abundance comes from SNIa, 
SNII, and SNIIx. The three different SNIIx models all 
produce similar results, and provide the silicon necessary to 
match the observational overabundance without unduly exac- 
erbating the low observed abundances of sulfur, argon and 
calcium. 

8.4. Discussion 

The anomalous abundance patterns found here in galaxy 
clusters are unique. The relative simplicity and uncomplicated 
nature of the X-ray emission from clusters and their role as a 
repository for all the enriched gas produced by supernovae 
makes them important objects for understanding the produc- 
tion and evolution of metals in the universe. 

However, no combination of SN la and SN II products us- 
ing current theoretical yields can produce the abundances 
observed with ASCA. Manipulating the standard models by 
changing the IMF or appealing to different physical models 
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Fig. 1 1. — A decomposition of the intermediate element abundances into 
their nucleosynthetic origins based on TablelTl The lowest band in each bar is 
the contribution from SN la, the middle band adds the yields from the canoni- 
cal SN II, and the top band includes the contribution from an early population 
of high mass, low metallicity Pop III supernovae progenitors. Iron and nickel 
both are dominated by SN la products, with no contribution by SN IIx from 
Pop III progenitors. 



such as delayed detonation can help mitigate only some of 
the inconsistencies between models and observations. A new 
source of metals is needed in order satisfactorily explain the 
cluster abundances. The low ratios of S/Si, Ar/Si and Ca/Si 
are indicative of a fundamentally different mode of heavy 
metal enrichment. We have identified three different SN mod- 
els that can produce these metals in the correct proportions 
and have calculated the number of events necessary to match 
the observations when combined with metals produced in or- 
dinary SNIa and SNII. These models all have in common 
very massive, metal poor progenitor stars. 

It is natural to associate these with the earliest generation 
of (Pop III) stars. This primordial population must have 
existed, and its distinct characteristics are well matched to 
those required to explain the ICM abundance anomalies: zero 
metallicity, and a top-heavy IMF that has an imp ortant con- 
tribut i on from very massive sta rs (Bromm. Copp i. & Larsonl 
119991 lAbel. Brvan. & NormanI |200 ) . Enhanced explosion 
energies are indicated, and nucleosynthetic production is high 
with an abundance pattern quite distinct from other SNII 
jHeger & Wooslev 2002) . 

nLoewensteinl (1200 ll) proposed hypernovae associated with 
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Pop III as a means of enhancing silicon relative to oxygen in 
the ICM, and inferred a hypernovae rate of the same order as 
the SNllx from Pop 111 stars derived here. The precise re- 
quired number is subject to assumed values of the ICM/stellar 
ratio, the IMF and SN 11 yields for Pop 11 stars, and the IMF 
and SNllx yields for Pop 111 stars, but is in the range of 
10-30 times larger than the average valu e (^Irrr ^4 x 10"^) 

Jredicted in the semi-analytic models of lOstriker & GnedinI 
1996). This could be the result of enhanced primordial star 
formation in these extremely overdense environments. Evi- 
dence also exists for more ne arly solar abundance ratios in 
less massive clusters (Finogue nov et al.l2 002). 

Both observations and theory suggest the idea that a large 
population of massive, metal-poor Pop 111 stars was present 
at very high redshift and constituted the first generation of 
stars. The heavy metal products from this first generation 
would have been formed at the same time as the first galax- 
ies or even before, and would now be widely dispersed in the 
ICM. Galaxy clusters serve as the largest repositories of both 
enriched gas from cluster galaxies and of gas expelled by the 
earliest generations of stars and then accreted onto clusters. A 
combination of X-ray observations well matched to observ- 
ing metal abundances in clusters and the importance of galaxy 
clusters as large retainers of Pop 111 enriched gas make these 
observations one of the best views onto the earliest genera- 
tions of stars in the universe. 

9. SUMMARY 

We have presented intermediate element abundances for 
galaxy clusters based on ASCA observations. Our measure- 
ments of the iron and silicon abundances agree w ith the 
past A SCA results of Fukazawal (Il997h and iMushotz kv et al.l 
il996h . but achieve much higher precision and extend the tem- 
perature range from 0.5-12 keV. The measurements of the in- 
dividual element abundances show some surprising new re- 
sults; silicon and sulfur do not track each other as a function 
of temperature in clusters, and argon and calcium have much 
lower abundances than expected. 

These results show that the a-elements do not behave ho- 
mogeneously as a single group. The unexpected abundance 
trends with temperature probably indicate that different en- 
richment mechanisms are important in clusters with different 
masses. The wide scatter in a-element abundances at a sin- 
gle temperature could indicate that SN models need some fine 
tuning of the individual element yields, or that a different pop- 
ulation of SN needs to be considered as important to metal 
enrichment in clusters. 

We have also attempted to use our measured abundances to 
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constrain the SN types that caused the metal enrichment in 
clusters. A first attempt to split the metal content into contri- 
butions from canonical SN la and SN 11 led to inconsistent re- 
sults with both the individual elemental abundances and with 
the abundance ratios of our most well measured elements. An 
investigation of different SN models also could not lead to a 
scenario consistent with the measured data, and we deduce 
that no combination of SN la and SN 11 fits the data. Another 
source of metals is needed. 

This extra source of metals must be able to produce enough 
silicon to match the measurements, but not so much sulfur, 
argon, or calcium to exceed them. An investigation of SN 
models in the literature led to three separate models that could 
fulfill these requirements. All three models were similar, and 
had as progenitor stars massive and/or metal poor stars. The 
combination of canonical SN la and SN 11 with one of the new 
models does a much better job of matching the observations. 
These sorts of massive, metal poor progenitor stars are exactly 
the stars that are supposed to make up the very early Popula- 
tion 111 stars. The conjunction of our required extra source of 
metals with SN from population Ill-like progenitors supports 
the idea that a significant amount of metal enrichment was 
from the very earliest stars. 

Clusters are a unique environment for elemental abundance 
measurements because they retain all the metals produced 
in them. The relatively uncomplicated physical environment 
in clusters also allows well-understood abundance measure- 
ments. Future abundance analyses using a large sample of 
XMM data will allow an even better understanding of the SN 
types and enrichment mechanisms important in galaxy clus- 
ters. 
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